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Abstract 

The systematic approach for the off-perturbative calculations in disordered systems is developed. 
The proposed scheme is applied for the random temperature and the random field ferromagnetic Ising 
models. It is shown that away from the critical point, in the paramagnetic phase of the random 
temperature model, and in the ferromagnetic phase of the random field one, the free energy contains 
non-analytic contributions which have the form of essential singularities. It is demonstrated that these 
contributions appear due to localized in space instanton-like excitations. 
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1 Introduction 

In very simplified terms, studies of classical statistical systems involves two main domains: (1) inves- 
tigation of the ground state, and (2) summation over fluctuations around this ground state. Although 
formally according to the definition of the partition function, one has to perform summation over the whole 
configurational space of a system, in reality it is never done. And it is not that we are doing something 
wrong. The point is that in most of the cases only very limited part of the configurational space which is 
relevant for observable thermodynamics. Very often, it the question, what this "relevant part" is (which 
involves the choice of the so called "relevant variables"), which is the most difficult. Studies of the sys- 
tems containing quenched disorder, in addition to the two items mentioned above, involves the third one 
(although, technically, very often it turns into the starting one), which is the averaging of self-averaging 
quantities over random parameters. Nevertheless, at a qualitative level, the situation here remains the 
same: only very limited part of the configurational space is relevant for observable thermodynamics. 

However, in some statistical systems, besides the ground state another local minimum (or minima) 
could exists. Let us consider extremely simplified situation, schematically shown in Fig.l, when in addition 
to the ground state the system has another local minimum located in the configurational space "far away" 
from the ground states, and separated from it by a big (compared to the temperature) energy barrier, 
which, however, remains finite in the thermodynamic limit. If we are dealing with the system which 
contains no quenched disorder, then the thermodynamic contribution due to this another state with an 
exponential accuracy will be simply of order of exp(— (3AE) (provided AE ^> T), where AE = E\ — Eq. 
The crucial point, however, is that to get the above exponential contribution, one has to know about 
existence of the other local minimum, otherwise, its contribution would be just missing. In other words, 
summing up the perturbation theory around the ground state, and even taking into account all non-linear 
terms of the Hamiltonian, responsible for the existence of the other local minimum, would not recover 
the contribution exp(— /3AE) of the other state, which is located "beyond barrier". It is these type of 
contributions which are usually called " off-perturbative" . 

In the studies of the effects produced by the quenched disorder, conditionally, one can distinguish 
two main domains of research: strongly disordered systems, like spin-glasses (where the disorder in the 
dominant factor), and the systems containing some kind of weak disorder which is supposed not to destroy 
the ground states properties of the corresponding pure system. Traditionally, magnetic statistical systems 
containing weak disorder, such as random bond ferromagnetic Ising models, are studied focusing mainly 
on modifications introduced into their critical properties at the phase transition point [1, 2]. In fact, as was 
pointed out by Griffith many years ago [3], modification of the critical behavior, is not the only qualitative 
physical phenomenon which can be produced here. 

Let us come back to the example shown in Fig.l. The presence of weak quenched disorder here, provided 
it does not ruin the global structure of the phase space, would just require supplementary averaging of 
the above exponential factor, exp(— (3AE) (since both the energy of the ground state Eq and the energy 
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Figure 1 : Schematic structure of the free energy landscape with two well separated thermodynamically 
relevant valleys 

of the excited state E\ are now the functions of the disorder parameters), but qualitatively, it would not 
modify the situation too much. Completely different and new physical phenomena comes into play when 
the structure of the phase space similar to that shown in Fig.l is created by the presence of randomness. 
In other words, this is the situation when weak quenched disorder, although it does not modify the ground 
state of the system, creates something completely new (absent in the corresponding pure system) , namely, 
the local minima states, somewhere at the periphery of the phase space, "far away" from the ground state 
of the system. 

According to the original observation by Griffiths [3] and later studies [4, 5], the presence of such off- 
perturbative states in the disordered ferromagnetic Ising model makes its free energy to be non-analytic 
function of the external magnetic field h in a whole temperature interval above the ferromagnetic phase 
transition point. Moreover, at least in some cases, this non-analyticity has the form of essential singularity 
in the limit h — > 0, [4,6]. It has to be stressed that all such contributions are just missing in the 
traditional (perturbative) RG treatment of the problem [2]. In more spectacular way the off-perturbative 
effects manifest themselves in the dynamical properties, producing the slowering down of the relaxation 
processes[7], as well as in the quantum systems (see e.g. [8] and references therein). 

Although at a qualitative level the origin of the off-perturbative contributions is more or less clear, their 
technical implementation, namely, the derivation of e.g. the non-analytic part of the free energy, turns 
out to be rather tricky problem. Usually, analytic calculations in disordered systems are performed in 
terms of the replica method. Many years ago Parisi has suggested [9], that the presence of additional local 
minima configurations in weakly disordered systems is related, in the replica approach, to the existence 
of localized in space and breaking replica symmetry instanton-like excitations (translation invariance and 
replica symmetry is recovered by taking into account all possible excitations of this kind). Later on 
there were several attempts of concrete realization of this idea for the random temperature [5] and the 
random field[10,ll] Ising models where it has been demonstrated that the corresponding saddle-point 
equations may indeed have instanton-like solutions. Next step has been done when the systematic method 
of summation over all such type of solutions, breaking symmetry in the replica vector order parameter, 
has been developed[12, 13]. In terms of this method the explicit form of the off-perturbative contributions 
in the random temperature Ising model has been derived[14]. 

In the present paper (following the recent study the original Griffith problem [6]) the systematic ap- 
proach for the off-perturbative calculations and its relation with the method of the vector replica symmetry 
breaking[12, 13] is formulated (section II). Then, proposed scheme is applied for the random temperature 
(section III) and the random field (section IV) ferromagnetic Ising models. It is shown that in both 
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systems at temperatures away from the critical point, namely in the paramagnetic phase of the random 
temperature model, and in the ferromagnetic phase of the random field one, the free energy contains 
non-analytic contributions which (as the functions of the disorder parameters) have the form of essential 
singularities. It is demonstrated that these contributions appear due to localized in space instanton-like 
configurations. Physical discussion of the obtained results is given in Section V. 



2 General scheme of calculations 

Let us consider a general (continuous) D-dimensional random system described by a Hamiltonian H [</>(x) ; £(x)] , 
where <j>{x) is a field which defines the microscopic state of the system, and £(x) are quenched random 
parameters. Let us suppose that in addition to the ground state, this system has another thermodynami- 
cally relevant (Griffith) region of the configurational space located "far away" from the ground state and 
separated from it by a finite barrier of the free energy (see Fig. I). In other words, it is supposed that the 
partition function (of a given sample) can be represented in the form of two separate contributions: 

Z = J V<p{^)e~ tm = e~ 0F ° + e~ 0Fl = Z + Z 1 (1) 

where F is the contribution coming from the vicinity of the ground state, and F\ is the contribution of 
the Griffiths region. Then, for the averaged over disorder total free energy we find: 

T = -±EZ = fb - [1 + Z 1 Z^ 1 ] (2) 

The second term in the above equation, which will be denoted by Fq, can be represented in the form of 
the series: 

Fg = —^ w - E (3) 



where 



Z n (m) = f[[ JJ f V+fl e-^-^.-.^.^.-.ej (4 ) 

6=1 ^ c=l 

is the replica partition function (H n [</>] is the corresponding replica Hamiltonian), in which the replica 
symmetry in the n-component vector field (f> a (a — 1, n) is assumed to be broken. Namely, it is supposed 
that the saddle-point equations 

SH n [0] 

^(xT = ' (a = 1 -< n) (5) 

have non-trivial solutions with the RSB structure 



</>i(x) for a = 1, 

/>:« = <! (6) 

^o(x) for a = m+ 1, ...,n 

with <f)i (x) 7^ 4> (x), so that the integration in the above partition function, eq.(4), goes over fluctuations 
in the vicinity of these components: 

= 0i(x) + <Pb(x), (6=1,..., m) 
^°>(x) - ^ (x)+xc(x), (c=l,...,n-m) (7) 

It should be stressed that to be thermodynamically relevant, the RSB saddle-point solutions, cq.(6), 
should satisfy the following tree crucial conditions: 

(1) the solutions should be localized in space, so that they are characterized by finite space sizes R(m); 
in this case the partition function, eq.(4), will be proportional to the entropy factor V/R D {m) (where V 
is the volume of the system), and the corresponding free energy contribution Fq, eq.(3), will be extensive 
quantity; 

(2) they should have finite energies E(m) = H n [<p*]; 

(3) the corresponding Hessian matrix of these solutions should have all eigenvalues positive. 
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Thus, in the systematic calculations one should find all saddle-point RSB solutions </>„(x), eq.(6), 
(satisfying the above three requirements), after that one has to compute their energies E(m) (for n — > 0), 
next one has to integrate over the fluctuations in the vicinity of these solutions, and finally one has to 
sum up the series 

Fg = ^ y tlT^ R -D {m) (p det f)~ 1/2 e-m** (8) 

p ^— ' m V /n=0 

where T is the (n x n) matrix 



<5 2 ff [0] 



d(p a d<p a > <t>=ct>* 



(9) 



Note that in the present approach the procedure of analytic continuation n — > is quite similar to that 
in the usual replica theory [15]: whenever the parameter n becomes an algebraic factor (and not the 
summation parameter, or the matrix size, etc.), it can safely be set to zero right away. 

The above scheme of calculations can be easily generalized for an arbitrary number of the Griffiths 
regions. For example, let us consider the situation which is qualitatively represented in Fig. 2, when in 
addition to the ground state, the system has two thermodynamically relevant Griffiths states. In this case 
instead of eq.(l) we will have 
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Figure 2: Schematic structure of the free energy landscape with three well separated thermodynamically 
relevant valleys 



Z = J Vc^{x)e- fm = e-< iFa + e-? Fl + e~ pF ' = Z + Z x + Z 2 (10) 
and correspondingly, instead of eq.(3) we find 

m—l 



' ' m=l fe=0 v ' 

~ lim y y m - Zn(Km) (ii) 

' m=l k=0 v ' 



Here, in the replica partition function 
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k r. m—kp n—m „ 

z n (k,m) = n / n / n / ^ (o) e-^^^^ ^ 

6=1 ^ c=l d=l 

the integration is supposed to be performed in the vicinity of the saddle-point replica vector 

(f>i (x), for o = l, fc 

0*(x) = < 02 (x), for a = k + 1, ...,m (13) 

</>o(x), for a = m + 1, ...,n 

(where </>i(x) 7^ </> 2 (x) 7^ 0o(x)) which is the solution of the saddle-point equations (5). Finally, for the 
Griffiths free energy contribution, instead of eq.(8) one obtain 

^— ^ ^ klim — k)\ V /n=o 

m=l r k=0 y ' 

where E(k, m) = i?„_>o [<p*] is the energy of a given solution, eq.(13), and T is the Hessian matrix, eq.(9). 



It is interesting to note that one can arrive to the same representations for the off-perturbative free 
energy contributions, eq.(8) or eq.(14), following the so called vector replica symmetry breaking scheme[12, 
13]. The starting point here is the standard replica definition for the averaged over disorder free energy, 

1 711 _ 1 

jr=_ i im (15) 

(3 rwO n 

where the replica partition function Z n = Z n is formally defined by the integration over all configurational 
space: 

n ,. 

z n = H v<j> a e-/» ff »[*i.-.*»] (16) 

0=1"' 

Now, let us suppose that in addition to the usual replica symmetric (RS) ground state configuration, the 
saddle-point equations (5) have another types of solutions, which are well separated in the configurational 
space from the RS state. In this case (again, denoting their contributions by the label "G") the replica 
partition function, eq.(16), can be decomposed into two parts: 

Z n = Z R s + Z G (17) 

Here Z^s contains all "routine" perturbative contributions in the vicinity of the ground state, and, as 
usual, (in the limit n — > 0) this partition function can be represented in the form: 

z RS = e-? nF ° (is) 

Thus, according to eq.(15) for the total free energy we get: 

T = F + F G (19) 

where 

F G = - lim -Lza (20) 

n->o fin 

contains all non-replica- symmetric contributions (if any). As an example, let us suppose that the saddle- 
point eqs.(5) have non-trivial solutions with the three groups structure like in eq.(13). Moreover, let 
us suppose that these solutions possess three crucial properties: (1) they are localized in space and 
characterized by finite spatial sizes R n (m, k); (2) they have finite energies E n (k, m); and (3) their Hessian 
matrices T a b, eq.(9), have all the eigenvalues positive. Than taking into account all possible permutations 
of the tree replica vector components the above free energy Fq can be represented in the form 



fG ,-ii4y u n - - y iu m[ - Z - (/MctTT'V^-) (21) 

n-o (in 4- ;m!(n-m)! k\(m - k)\ R%(k,m) V J 
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To perform the analytic continuation n — > in the above expression the parameter n must enter as an 
algebraic factor, and not as the parameter of summation. This can be achieved if we represent the above 
series in the following way: 

F ° - ~ iS Tn £ Tin, + l)'r"(I -m + l) £ M^- »)■ R S l, m ) (*» *)"" ^ 

Here the summation over m is extended beyond m = n limit since the gamma function T(z)is equal to 
infinity both at z — and at all negative integers. Now using the relation: 

r(-«) = , ^ , ; (23) 

V ; zT(z) sin(7rz) V ' 

and referring to "good" analytical properties of the Gamma functions, we can perform the analytic con- 
tinuation n — > 0: 

T(n+1) T(n+1) 



r(m + l)r(n - m + 1) r(m + l)r[-(m - 1 - n)] 

r(n + l)(m — 1 — n)r(m — 1 — n) sin[7r(m — 1) — 7rn] 
7iT(m + 1) 



sin(Trn) ^_^ m -i T (n + 1)T (m - n) 



7T v ' r(m + l) 



~ ^ (24) 

m 

Substituting this into eq.(22) we obtain eq.(14) (where i?(fc, m) = R n =o(k, m) and -E(fc, m) = E n= o(k, m)) 
Now, in the next two sections I am going to demonstrate how the above general scheme works in the 
concrete cases of the random temperature and the random field ferromagnetic Ising models. 

3 Random temperature Ising model 

Let us consider weakly disordered D-dimensional Ising model described by the continuous Ginzburg- 
Landau Hamiltonian: 



H 



\ (V^(x)) 2 + \ (r - 5r(x)) </> 2 (x) + i# 4 (x) 



(25) 



The disorder is modeled here by a random function r(x) which is described by the Gaussian distribution, 

P[St] ^ exp^-i- J d D x(5T(x)) 2 ^ , (26) 

where u is the small parameter which describes the strength of the disorder, and po is the normalization 
constant. We are going to consider this system in the paramagnetic phase away from the critical point, 
so that the reduced temperature parameter r will be taken to be positive and not too small (it will be 
demonstrated below that in dimensions D < 4 the limitation on the value of r is given by the usual 
Ginzburg-Landau condition, r 3> g 2 ^ 4 ^ D ^). 

Weakly disordered systems described in terms of the continuous Ginzburg-Landau Hamiltonian have 
been usually studied in the framework of the renormalization-group (RG), (perturbative) treatment. In 
this approach one is able to perform the systematic integration over all fluctuations at the background of 
the homogeneous state up to the scales of the correlation length R c (r). In some cases, this makes possible 
to derive the leading singularities of the thermodynamical functions in the critical point, at r — > 0, where 
R c (t) diverges [2]. 

However, considering paramagnetic phase of this system, intuitively, it is clear that the contributions 
of non-homogeneous local minima configurations at scales bigger that the correlation length, which exist 
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due to rare localized in space " ferromagnetic islands" with negative effective value of the local temperature 
parameter (r — St) are missing in the traditional RG treatment of the problem. To what extend these 
states are relevant for the critical behavior (at r — > 0) is still unclear[16]. In the present study , however, I 
am going to address much more simple question: what is the explicit form of the free energy contributions 
due to such off-perturbative states away from T c , (at t>t 9 ). 

In fact, at purely heuristic level, it is not so difficult to estimate form of these contributions. Let us 
consider the spatial island of the linear size L characterized by the typical value of the " local temperature" 
(r — St) = — £ < Its probability is exponentially small, 



(27) 



and therefore such islands arc well separated from each other and can be considered as non-interacting. 

It has to be noted that the island with small (negative) value of the local temperature parameter 
— £ can be characterized by the mean-field "up" and "down" states only if its size is bigger than its 
local correlation length i? c (£) ~ £ -1 ^ 2 - Thus, the contribution to the free energy coming from the local 
ferromagnetic states of such islands with the exponential accuracy can be estimated by their probability: 



F G 



d£ / dLexp 



d£exp 



- {const) tL + tl r D/2 



(28) 



Here in the integration over £ the leading contribution comes from the vicinity of the saddle-point value 

D 



A-D 



(29) 



(which is positive in dimensions D < 4, and 3> r g provided t 3> T g ). In this way we obtain the following 
estimate for the off-perturbative contributions coming from rare locally ferromagnetic islands: 



F G 



exp 



- (const) - 



r (4-D)/2 



(30) 



Now let us consider how this result (including the value of the (const) factor) can be derived analytically 
in terms of the systematic approach developed in the previous section. This derivation has been already 
reported elsewhere [6, 14]. Here I am going to give some more details about the corresponding replica 
instanton solutions, but in general this section can be considered just as a "warming up" exercise before 
passing to more difficult calculations for the random field model considered in the next section. 

Performing the standard Gaussian integration over random parameters <5r(x), for the replica partition 
function one gets 



z n = f[ Jvm*) e~ 0HnW 



where 



H n [cj>] 



a=l 



Z a=l 



1 

a— 1 



1 n . n 



a=l 



a,b=l 



(31) 



(32) 



is the corresponding replica Hamiltonian. The saddle-point configurations of the fields <f> a (x) are defined 
by the equations 



-A<f> a (x) + T(t> a (-x) + 3 <^(x) - u0 a (x) 



>g(x) 







(33) 



Below we are going to demonstrate that besides the trivial solution o (x) = these equations have 
non-trivial localized in space instanton-likc solutions with the RSB (two groups) structure: 




for a - 



, m 



for a = m + 1, n 



(34) 
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Substituting this anzats into the saddle-point eqs(33) and into the Hamiltonian, eq.(32), we find that (in 
the limit n — > 0) the instanton configuration <pi {x) is defined by the equation 



- A<£i (x) + rfa (x) - A(m)0i (x) 3 = 
which is controlled by the parameter 



A(m) = um 



and the energy of this configuration is 

E{m) = m J d D x 



; (V« 2 + ^->(m)^ 



(35) 
(36) 

(37) 



In what follows the parameter A(m) will be assumed to be positive. In other words, the solution, which 
we are going to derived below, exists only for m such that m > [g/u] (where [...] denotes the integer 
part). It has to be noted that one should not be confused by the "wrong" sign of the coupling 4 term 
in the above equations. In fact, it can be shown that the integration over the replica fluctuations around 
considered solution in the limit n — > yields the Hessian matrix which has all the eigenvalues positive 
(this is the usual situation for the replica theory, where the minima of the physical quantities in the limit 
n — > turns into maxima of the corresponding replica quantities[13, 15]). 
Rescaling the fields, 

01 (x) = 



A(ra) 



(38) 



and introducing z = x-/r, instead of eq.(35) one get the differential equation which contains no parameters: 

-Aip(z) +V(z) - ^ 3 (z) = (39) 
Correspondingly, for the energy of this configuration, eq.(37), one obtains: 



- T {4-°)/*Eo(D) 



E{m) = -JUL 
um — g 

where the quantity E (D) depends only on the dimensionality of the system: 



E (D) = J 



l(V^(z)) 2 + ^ 2 (z)-^ 4 (z) 



(40) 



(41) 



It can be shown (see e.g. [17]) that in dimensions D < 4 eq.(39) has the smooth (with V'(0) = 0) 
spherically symmetric instanton-like solution ij)(r) (where r = |z|) such that: 



ip{r < 1) 



^(r = 0) = Vo ~ 1, 



ip{r > 1) 



0. 



(42) 



The energy E (D) of this solution is a finite and positive number. In particular, in dimensions D = 3, 
ip ~ 4.34 and E Q ~ 18.9 (see Fig.3). 



As the dimension parameter D approaches the upper critical dimensionality D c = 4, from below the 
value of the field at the origin tpo(D) tend to infinity (see Fig. 4), while the energies of the corresponding 
instanton configurations Eq{D) approach the finite universal value E^(D — > 4) = E* ~ 26.3. Above 
dimensions D = 4 eq.(39) has no smooth instanton-like solutions. In other words, described in terms of 
the dimensions parameter D, when passing the critical value D c = 4 from below, the instanton solution 
disappears in the discontinuous way. 
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Figure 3: The instanton solution ip(r) of eq.(39) in dimensions D = 3 



Note that according to the rescaling, eq.(38), the size of these instanton solutions in terms of the 
original fields <j>i(x) is R c (t) = t -1 / 2 (which is the usual correlation length of the Ginsburg-Landau 
theory) and it does not depends on m. Note also that due to obvious symmetry property 4> a — ► — 4> a of 
the original saddle-point eqs.(33) (which is valid for all non-zero replica field components independently), 
the above instanton solution </>„(x), eq.(34), has additional degeneracy factor 2 m . 

The final step is the integration over fluctuations <£> a ( x ) at the background of the above instanton 
solution. Substituting a (x) = 0*(x) + i£ a (x), and expanding the Hamiltonian up to the second order 
in <p a (x), one has to perform the standard Gaussian integration. These calculations, although slightly 
cumbersome, are quite straightforward (for the details see Refs.[6, 14]). In the result for the Hessian 
factors one gets 



(detTj 



-1/2 

~ exp 

rwO 



3m , 9 

Mo 



2(wm — g)' 



(43) 



Comparing this with factor exp[— E(m)], where E(m) is the instanton energy, eq.(40), we see, that under 
condition 

r » r g = g 2 '^ (44) 

the contribution of fluctuations can be neglected. This is not surprising because eq.(44) is nothing else, 
but the Ginzburg-Landau criterion which defines the temperature region away from T c , where the critical 
fluctuations are irrelevant. On the other hand, it has to be stressed that in the close vicinity of T c (at 
T < T g)i where the critical fluctuations are relevant, the Gaussian approximation used for obtaining the 
result, eq.(43), can not be valid anymore, and to derive the corresponding fluctuations contribution one 
would have to start some kind of RG procedure which would properly take into account non-Gaussian 
interactions. Thus, the above result for the fluctuations contribution, eq.(43), either can be considered as 
the small correction (at r 3> r s ), or otherwise (if it is not small), it is not valid (at r < r g ). 

Thus, substituting the value of the instanton energy, eq.(40), its size R = r^ 1 / 2 as well as its degeneracy 
factor 2 m into the series, eq.(8), we get 



F G ~ -Vt d ' 2 Y LJ2 2 m exp 

j—f-, m 

m=[g/u] + l 



. Eo{D) ^_ T ^-D)/2 

um — g 



(45) 



The exact summation of this series seems to be rather tricky problem, but with the exponential accuracy 
it can be estimated in a very simple way. One can easily see that in the limit of weak disorder, at u <C g, 
the leading contribution in this summation comes from the region m 3> g/u^> 1 (where the exponential 
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Figure 4: The instanton solutions ip(r) of eq.(39) in various dimensions: (a) D = 2, ipo — 2.21; (b) D = 3, 
ip ~ 4.34; (c) D = 3.5, Vo ^ 7.92; (d) D = 3.8, Vo ^ 16.4; (e) D = 3.9, ip ^ 26.4; 



factor in eq.(45) becomes m-independent) and this contribution is 

/ (4-£>)/2\ 

F G ~ exp f -£ (I>) J (46) 

We see that obtained off-perturbative (Griffith- like) part of the free energy, as the function of the disorder 
parameter in the limit u — > 0, has the form of the essential singularity. Note again that this contribution 
exists only in dimensions D < 4. As discussed above, at D — > 4 (the upper critical dimensions), the 
dimensionless instanton energy factor Eq(D) approaches the finite universal limiting value E„ ~ 26.3. In 
our world, in three dimensions, E (D — 3) ~ 18.9. 



4 Random field Ising model 



To study the off-perturbative effects in the random field Ising model we are going to use again the 
Ginzburg-Landau continuous representation: 



H 



d u x 



\ (V0(x)) 2 



ir0 2 ( X ) + i# 4 (x)-^)0(x) 



Here the random function h(x) is described by the Gaussian distribution, 



P[/i(x)] =p exp 




d D xh 2 (x) 



(47) 



(48) 



where ho is the small parameter which describes the effective strength of the random field, and po is 
the normalization constant. Unlike the random temperature model, considered in the previous section, 
here we are going to consider the system in the low-temperature ferromagnetic phase (supposing that the 
dimensionality D is such that this phase exists), so that the reduced temperature parameter r will be 
taken to be negative, r = — |r|. Again, we will place the system away from the critical point, assuming 
that the absolute value |r| is not too small. As usual, to avoid the effects of the critical fluctuations (in 
dimensions D < 4) we impose the condition |r| 3> g 2 ^ 4 ~ D \ 
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Let us suppose that in the absence of the random fields the ferromagnetic ground state of the system, 
eq.(47), is "up". This state (at the mean- field level) is characterized by the order parameter 




and the energy density 



eo 



Mi 



(49) 



(50) 



In the usual perturbative approach the effects produced by the random field term of the Hamiltonian, 
eq.(47), together with the thermal fluctuations could be calculated in the systematic way in terms of the 
RG procedure (see e.g. [18] and references therein). This approach is designed to take into account all 
degrees of freedom at scales less that the correlation length, R c (t) ~ |r|~ 1//2 . 




Figure 5: Two alternative local minima field configurations in a spatial "island" of the linear size L, where 
the average value of the random field h is negative and its absolute value is not too small. 



On the other hand, at scales bigger that the correlation length we can observe completely different 
type of thermal excitations. Let us consider a spatial island of the linear size L, where the average value 
of the field h is negative and its absolute value is not too small. Then, in addition to the state "up" (with 
slightly modified value of the order parameter), another local minimum with orientation "down" can exist 
in this island (Fig. 5). To be stable, the gain in the energy due to the interaction with the field, 

E h ~ -L D \h\4> (51) 

should overrun the loss of energy due to the creation of the domain wall, 



L 



D-l<t>l 



R, 



Thus, such double-state situation in the considered island is created provided 

*o Irl 



\h\ > 



= h r 



(52) 



(53) 



LR C LJg 

According to eq.(48) the probability to find an island of the size L with the average value of the field h is 

h 2 



P(L, h) ~ exp 



2h? 



L 



(54) 



Then the contribution to the free energy of such rare "flipped" states can be estimated by their probability: 

h 2 



F G 



dL / dh exp 



- D 



2h 2 



dL exp 



R c 



<amSt) k LD ~ 2 . 



exp 



- {const) - 



h l9 



(55) 
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Note that to obtain this result in the above integration over L the dimensionality of the system D must 
be bigger than two (otherwise the integral will become divergent). This is nothing else but slightly 
modified version of the good old Imri-Ma arguments[19] which tells that in dimension D < 2 flipping of 
magnetizations in big spatial islands can become energetically favorable, which indicate the instability of 
the global ferromagnetic state. Here we assume that the ferromagnetic state is stable, and we see that 
rare off-perturbative flipping excitations produce non-analytic contribution to the free energy, which in 
the limit ho — > has the form of essential singularity. 



Now we are going to re-derive the above prediction, eq.(55), in terms of much more rigorous systematic 
procedure described in section II. Coming back to the original Hamiltonian, eq.(47), after the Gaussian 
averaging of the replicated partition function over the random function /i(x), one obtains the replica 
Hamiltonian 



d° x \ E (v^) 2 -\\r\Y.€ + \ 9 E - \ h l E 

. a=l o=l o=l a, 6=1 

The saddle-point configurations of the fields (f) a {x) are defined by the equations 

-A0 o (x)-|r|0 o (x)+ fl ^(*)-fco (X>( X )J =° 



(56) 



(57) 



vb=l 



Below we are going to demonstrate that besides the obvious (replica symmetric) ferromagnetic solution 
<j> a (x) = y/\r\/g these equations have non-trivial localized in space instanton-like solutions with the RSB 
two-groups structure: 



C(x) 



^(x-y/jrl) fora=l,...,m 
^o(x-\/]t[) for a = m + 1, . 



(58) 



Substituting these rescaled fields into the saddle-point eqs(57) and into the Hamiltonian, eq.(56), we find 
that (in the limit n — > 0) the instanton configuration {V>i(z), ipo( z )} (where z = x-y/jrf) is defined by the 
two equations 



-Aip! -ipi+ipf- A(m) (ipi - ip ) = 



-AV>o - Vo + - x (m) (ipi -ipo) = (59) 



and its energy is 

\ T \2-D/2 r 

E{m) = m^- J 

where 



d D 



\ [(v^o 2 - (v^o) 2 ] -\Wi- ri] + V'o 4 ] - \Km) [Vi - M 2 



(60) 



A( TO ) = (61) 



We are looking for the localized in space (spherically symmetric) solutions of the eqs.(59), such that the 
two functions ip\{r) and V'o( r ) (where r = |z|) are different from each other in a finite region of space, and 
at large distances they both sufficiently quickly approach the same value ip = 1, so that the integral in 
cq. (60) will be converging. Simple analysis of the structure of the " potential energy" 

u(A,i>o) = +\[1>\-<] ~\^i-M 2 (62) 

shows that until the parameter A is small (so that the last coupling term in the above expression 
is just a small correction), the potential U(ip\,ipo) has 9 saddle-points (in the vicinity of the points 
(0; 0), (0; ±1), (±1; 0), (±1, ±1) and (±1;=f1)). In this situation the two fields ip\ and ipo are effectively 
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independent, and no instanton-like solutions described above can exist. When increasing the parameter 
A, starting from 

A > A c ~ 0.23 (63) 

only 5 saddle points of the potential /7(V>i, V'o) remains in the plain (ipi;ipo). They have coordinates: 
(0;0),(±1;±1) and (±ipl(X);±ip^(X)), where < ipl >0 < 1 (in particular, ipl(X c ) ~ 0.17 and %{X C ) ~ 
0.90). It is crucial that at the points (±^>* ; ±^5) * ne potential U(ipi,ipo) has the maxima. It is due to the 
existence of these maxima that at A > A c the instanton solutions become possible. 
Let us consider the limit A(m) S> 1, or 



C hl 



+ 1 



(64) 



In this limit, according to eqs.(59), the two fields ipi and tpo must be close to each other. Redefining, 

1 



V>i(r) 



(65) 



in the leading order in A 1 instead of eqs.(59) we get much more simple equations: 



-Aip - ip + V> - 2x 



-A x +(3^ 2 -l)x = (66) 

which contain no parameters. For the energy of the configurations described by the two fields ip{r) and 
X(r) instead of eq.(60) (again, in the leading order in A -1 ) we find the value, which does not depend on 
the summation parameter m, 

E = ^ E ° (D) (67) 

where 

E (D) = J d D z [(W)(Vx) + (^ 3 - V)x - X 2 ] (68) 

is the universal quantity which depends only on the dimensionality of the system. 

It turns out that in dimensions D < 3, the system of eqs.(66), indeed has smooth instanton-like 
spherically symmetric solution which has finite and positive energy E (D). Within the limited spatial 
region r < r c ~ 1, the values of the fields ip{r) and x(r) are finite and of the order of their values at 
the origin, ipo ~ 1 an d xo ~ 1- O n the other hand, at r > r c the function exponentially quickly 
approaches 1, while the function \{ r ) exponentially tents to zero. The illustration of the instanton solution 
in the dimension D = 2.9 is given in Fig. 6, where ipo — —0.818, xo — —0.284 



Fig. 7 demonstrates the corresponding "trajectories" of the instanton solutions in the plane (ip,x) a t 
various values of the dimension. As the dimension D approaches the value D c = 3 from below, the starting 
values ^o^—l an d Xo — ► 0. Above three dimension the instanton solution disappears. Thus we have to 
conclude that D = 3 is the upper critical dimension for the considered Griffiths phenomena in the RFIM, 
in agreement with the earlier suggestion [10, 11] as well as with the recent studies of similar instanton-like 
configurations in the presence of external magnetic field [20] 



Let us come back to the general expression for the off-perturbative part of the free energy, the series, 
eq.(8), where the summation over m starts now from m = m c , eq.(64). Noting that the instanton energy 
E(m) is the decreasing function of m, we can conclude that with the exponential accuracy this converging 
series can be estimated by its asymptotic part at ra » m c . Thus, substituting here the value of the 
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Figure 6: The instanton solution ip(r), x( r ) °f eqs.(66) in dimensions D = 2.9 



instanton energy, eq.(67), (and neglecting the critical fluctuations), with the exponential accuracy we get 
the result 



which perfectly agrees with the " hand- waving" estimate, eq.(55). This non-analytic in ho contribution, 
which has the form of essential singularity at ho — > 0, is valid only in dimensions D < 3, and at temper- 
atures not too close to the critical point, at |t| ^> g 2 ^ 4 ~ D K The dimension D = 3 is marginal for this 
kind of phenomena. Therefore the investigation of the Griffith-like contributions in the three-dimensional 
RFIM requires much more discreet analysis. 

5 Discussion 

In this paper the systematic method for computing off-pcrturbative thermodynamic contributions in disor- 
dered systems has been proposed. It has been tested on two the most popular classical statistical systems 
containing quenched disorder: the random temperature and the random field Ising models. In both cases 
the off-perturbative contributions as the functions of the parameters which describe the strength of the 
disorder have the form of the essential singularities, eqs.(46), (69). 

Of course, thinking about possible experimental or numerical tests, the validity of the obtained results 
is rather limited. On one hand, since the consideration has been done in terms of the continuous Ginsburg- 
Landau Hamiltonian, one has to place the system sufficiently close to the phase transition point, so that 
the correlation length would be large compared to the lattice spacing. On the other hand, since the 
present study completely neglects the critical fluctuations, the system has to be sufficiently far from the 
critical point. Formally, in terms of the Ginsburg-Landau Hamiltonian, these two requirements can be 
easily satisfied: it is sufficient to demand that (1) the coupling parameter g is small, and (2) the reduced 
temperature parameter r is bounded by the condition g 2 /^ D ) -C |r| <C 1 (where D is the system 
dimensionality, D < 4). Besides, the strength of the disorder must be small: u <C g in the random 
temperature model, and h$ <C vM m ^ nc ran dom field one. 

On the other hand, keeping in mind more general perspectives, the proposed approach, in my view, 
may open the way to study the nature of the phase transitions in the considered systems. It is generally 
believed that it is the off-perturbative states which makes the study of the phase transition in the random 
field Ising model so difficult [18]. It is remarkable, that according to the present study, quite similar off- 
perturbative contributions are also present in the random temperature Ising systems where, at least at 
the qualitative level, the nature of the phase transition was traditionally believed to be well understood 



F G ~ exp -E (D) 




(69) 
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Figure 7: The instanton solution "trajectories" in the plane (ip,x) m various dimensions: 
(a) D = 2, i[>o ~ 0.001, xo - -0.5173; (b) L> = 2.5, - -0.325, X o - -0.5776; (c) D = 2.7, 
?/>o ^ -0.520, xo - -0.5300; (d) D = 2.9, - -0.818, xo - -0.2844; (e) D = 2.92, - 
-0.988, xo - -0.0212. 



(see e.g. [2, 13]). This supports recent suspects[16] that the off-perturbative effects could be quite relevant 
for the critical properties of the disordered ferromagnetic systems. 

Another interesting observation is that in terms of the considered off-perturbative contributions, the 
situation, when approaching T c from above and from below, looks totally asymmetric both in the random 
field and in the random temperature systems. The considered effects are present in the ferromagnetic 
phase, at T < T c , of the random field model, while they are absent in its paramagnetic phase at T > T c . 
On the other hand, the situation in the random temperature model is similar, although "reversed": the 
off-perturbative contributions are present in the paramagnetic phase, at T > T c , while they are absent in 
the ferromagnetic phase, at T < T c , May be it is this asymmetry, which makes the nature of the phase 
transitions in these systems to be so non-trivial? 

One more qualitative observation is that according to the present study, the off-perturbative contribu- 
tions are absent in the random temperature model at dimensions D > 4, and in the random field model 
at dimensions D > 3. As for the random temperature systems, this is not surprising: all the previous 
studies were definite that at D > 4 the disorder is irrelevant for the phase transition. What does this 
mean for the random field model is much less clear, because here it is well established that its upper 
critical dimensionality is equal to 6 (the dimensionality above which the critical behavior is described by 
the Gaussian theory, and the presence of the random fields is irrelevant for the phase transition). Well, 
of course, the absence of the off-perturbative contributions in dimensions 3 < D < 6 does not mean, that 
the random fields are irrelevant. Probably it indicates that here all the random field effects can be taken 
into account in the framework of the perturbative RG procedures. 

To answer all the above questions (as well as many other important questions which were not formulated 
here), the only thing which remains to be done is to find the way to overcome the Ginzburg-Landau 
limitation |t| ^> g 2 /( 4 ~ D \ and to take the limit t — > 0. To do that one has just to formulate a theory 
which would properly take into account the critical fluctuations on top of the instanton-likc configurations 
described in this paper. 
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